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Abstract 



The anomalous acceleration ap of the Pioneer 10/11 - the Pioneer effect - 
has remained unexplained. We suggest an explanation based on the interac- 
tion of the spacecraft with a long-range scalar field, (/>. The scalar field under 
consideration is external to gravity, coupled to the ordinary matter and under- 
Q\ . goes obedience to the equivalence principle. In addition to its self-interaction 

Q I term, the field is determined by an external source term proportional to the 

O^' Newtonian potential in the weak fields limit it result a long-range accelera- 

5-H . tion ap, asymptotically constant within the region of the solar system hitherto 

^^' crossed by the spacecraft. Also, the limit 0.1 10^^ cm/s^, which follows from 

the Viking ranging data, is satisfied in the region of the terrestrial planets 
S^ . (in particular at the positions of the Earth and Mars), for a (/)-field of mass 

H I 771,0 > 1.8 10~^^ eV/c^. The proposed solution gives the correct order of 

magnitude for ap, as observed so far, and predicts the decline of ap in the 
form of damped oscillations beyond 97 AU. An estimate of the cosmological 
constant is also made by taking into account the contribution of the vacuum 
energy density of the scalar field in galactic dynamics and particularly in the 
outskirts of the dwarf galaxy DDO 154. 



1 Introduction 

Recently, results from an almost twenty years study of radio metric data from Pi- 
oneer 10/11, Galileo and Ulysses spacecraft have been published by a team of the 



NASA (Anderson et al. |1|), indicating an apparent anomalous, constant, accelera- 
tion acting on the spacecraft with a magnitude of the order 8.5 10~® cm/s^, directed 
towards the Sun, to within the accuracy of the Pioneers' antennas and a steady fre- 
quency drift, a "clock acceleration", of about —6 10^^ Hz/s. A number of potential 
causes have been ruled out by the authors, namely gravity from Kuiper belt, gravity 
from the Galaxy, spacecraft "gas leaks", anisotropic heat (coming from the RTGs) 
reflection off of the back of the spacecraft high-gain antennae (Katz's proposal [0, 
see Anderson et al. 0|), radiation of the power of the main-bus electrical systems 
from the rear of the craft (Murphy's proposal 0], see Anderson et al. 0), errors in 
the planetary ephemeris, and errors in the accepted values of the Earth's orientation, 
precession, and nutation, as well as nongravitational effects such as solar radiation 
pressure, precessional attitude-control maneuvers and a possible nonisotropic ther- 
mal radiation due to the Pu^'^^ radioactive thermal generators. Indeed, according 
to the authors, none of these effects explain the apparent acceleration and some 
are 3 orders of magnitude or more too small, so they conclude that there is an 
unmodeled acceleration towards the Sun of (8.09 ± 0.20) 10~^ cm/s^ for Pioneer 10, 
(8.56±0. 15)10"® cm/s^ for Pioneer 11, (12±3)10"® cm/s^ for Ulysses and (8±3)10-® 
cm/s^ for Galileo. The authors plan to utilize two different transmission frequencies 
in further analysis to give an answer to whether there is some unknown interaction 
of the radio signals with the solar wind. 

Since no "standard physics" plausible explanations for the residual acceleration 
has been found so far, the authors considered the possibility that the origin of the 
anomalous signal is the effect of a modification of gravity, for instance by adding 
a Yukawa force to the Newtonian or Milgrom's proposed modification of gravity 
(Milgrom 0). They concluded however that neither easily works. 

If the cause is dark matter, the amount needed to be consistent with the accuracy 
of the ephemeris should be only of order a few times 1O~^M0 even within the orbit 



of Uranus (Anderson et al. |0). Above all, the authors point out that the residual 
acceleration is too large to have remained undetected in the planetary orbits of 
the Earth and Mars. Indeed, the Viking ranging data limit any unmodeled radial 
acceleration acting on the Earth and Mars to no more than 0.1 10^^ cm/s^. Because 
of this severe constraint, the authors argue that, if the anomalous radial acceleration 
is of gravitational origin, it probably violates the principle of equivalence. But, an 
alternative is that the anomalous acceleration is asymptotically constant, rather 
than constant at all radii from the center of the solar system to the present location 
of Pioneer 10/11 spacecraft. 

In this paper we propose an alternative explanation, based on the possible exis- 
tence of a long range (non gravitational) scalar field, 0, which respects the (weak) 
equivalence principle. This possibility was previously introduced by Mbelek p|, to 
account for the rotational curves of spiral galaxies, as an alternative for dark matter. 
It gives, for the Pioneer 10/11 spacecraft, the correct order of magnitude for both 
the anomalous acceleration, ap, and the clock acceleration, at. It proves to remain 
consistent with the planetary orbits determined from the Viking data. 

As for the ordinary matter, the 0-field is a gravitational source through its 
energy-momentum tensor. A forthcoming paper will present the fundamental sym- 
metry that may support it, from the background of classical fields theory. The plan 
of this paper is as follows : in section 2, we set the 0-field equation. Then, after 
linearization we divide space in three characteristic regions and find approximate 
exterior solutions for a static spherically symmetric source, actually the Sun. In 
section 3, Einstein equations are solved in the weak fields approximation to account 
for the metric tensor in the presence of the 0-field (out of the Sun). In section 4, 
the equation of motion of a test body in the presence of the 0-field is established. 
Solutions are found in the weak fields and low velocity limit. Then the anomalous 
long-range acceleration ap is derived for the different regions of space. In section 5, 



an interpretation of the data is proposed. In section 6, the steady frequency drift at 
is derived by using the equivalence principle. We finally conclude by an estimation 
of the cosmological constant by exploiting the declining part of the rotational curve 
(RC) of the dwarf galaxy DDO 154. 

2 The scalar field equation 

A manner to generate a constant radial acceleration could result from the introducion 
of a linear potential term in the Lagrangian of a test particle. An example is provided 
by the exterior solution of the locally conformal invariant Weyl gravity for a static, 
spherically symmetric source (Mannheim and Kazanas [^). Unfortunately, Per lick 
and Xu [|1^], by matching the exterior solution to an interior one that satisfies the 
weak energy condition and a regularity condition at the center, show that this leads 
to contradiction of Mannheim and Kazanas's suggestion. They conclude that the 
conformal Weyl gravity is not able to give a viable model of the solar system. 
This paper presents an alternative solution, under the form of a real scalar field, 
external to gravity but which satisfies the equivalence principle. We show below that 
it leads to the desired " Pioneer effect ", although it does not modify, as required, 
the orbital properties of the internal planets. The field obeys the equation 

V.V> = -t/'(0) - J, (1) 

where the symbol V^ stands for the covariant derivative compatible with the Levi- 
Civita connection. Equation (|I|) may be derived from Einstein equations provided 
that the energy-momentum tensor of the 0-field is of the form T^^^ = d^(f) dycj) — 
g^lv[\^\(t) d^(f)—U{(f>)~J Jd(j)] (up to a positive multiplicative dimensionality constant, 
K, for is dimensionless in this paper) and the energy-momentum tensor of the 
ordinary matter (matter or radiation other than the 0- field) is divergenceless (e.g., 
zero for the exterior solution and of the perfect fluid form for the interior solution). 



In the rest of the paper, we apply the weak field approximation to the real classical 
scalar field and to the gravitational potentials, so that Newtonian physics apply. 
For a weak gravitational field, equ(P above will write merely 

d,d^ct> = -U'{<P) - J. (2) 

The potential U denotes the self-interaction of (j), and we note U'{(j)) = ^. The 
source term J, an external source function, takes gravity into account, as a source 
for the field (p. Of course, also acts as a source for gravity (through Einstein 
equations). We consider this latter action below in section 3. 

In the weak field approximation, J depends on the Newtonian gravitational 
potential V^, the only relevant scalar quantity related to a weak gravitational field. 
Thus, we write at first order J = J(-^) ~ — Y^, where the constant tq defines a 
characteristic length scale (see subsection 5.1.1 for an estimation of ro). The minus 
sign comes from the requirement that the effect of is similar to that of gravitation, 
so that -£: and ^^ have the same sign (as we will see, the 0-field generates an 
acceleration term ^, up to a positive multiplicative factor). This is in accordance 
with our previous study on the RC of spiral galaxies in which we found the positivity 
of -£: necessary for the 0-field mimics a great part of the missing mass |^. In the 
solar system, -^ remains positive. The scalar field is positive definite throughout 
this paper. 

Here we explore the effect of in the solar system, i. e., in the potential Vn = 
—c^ rs/2r created by the static central mass of the the Sun, r being the radius from 
the centre (we choose as usual a zero value of the Newtonian potential at infinity), 
and Ts the Schwarzschild radius of the Sun. The problem has spherical symmetry, 
so that equation (H) yields finally 

dr^ r dr 2r Tq 



We will calculate the resulting 0-field and show that it creates an asymptotically 
constant acceleration: we solve the equation with the limiting condition (imposed 
by the weak fields approximation) that the field and its derivative -£: are bounded 
for any given region of space. In addition, the field (p must vanish (up to an additive 
constant) if one sets M = 0, since the central mass M is its source ; the same 
condition applies to -£;. 

As a first step to resolve equ(^, let us neglect for the moment the contribution 
of the self-interaction. The solution is 



r^ A 



C + -^r - - (4) 



and thence, 



— = — + — (5) 

dr Arl r^' 



where A and C are constants of integration. The constant of integration A of the 
dimension of a length obviously depends on r^ since the central mass is the source of 
the 0-field. Accordingly, we may set A = C'^s/2, where ^ is a positive dimensionless 
constant that we will assume hereafter of the order unity. The positivity of ( is 
inferred by the positivity of the spatial derivative -£: at any distance from the centre. 
Note that the particular value C = 1 involves the identity of the potential term —A/r 
with the Newtonian one V^/c^. We will show that, in a certain radius range and 
at sufficiently large distances from the centre, this represents the true solution to 
equ(0): the radial acceleration induced by the (/)-field (neglecting its self-interaction), 
proportional to -£:, remains asymptotically constant. This will be referred to as the 
"Pioneer effect" throughout. 

In order to solve the complete equation, we need to know the form of f/(</>), although 
we will see that many results remain independent of this choice. To illustrate, we 
choose here a quartic self-interaction potential U = f/(0) + |/i^0^ + |0^, where 
0" < is the self-coupling coefficient of the scalar field ; ^/i^^^ is the "mechanical" 



mass term with fi = ^^^, m^ denoting the mass of the scalar field. The reason for 
choosing a quartic polynomial form is that the corresponding quantum field theory 
should be renormalizable (Madore p!I|). 

This potential presents two extrema : one minimum at 0/ = 0(r/), with U'{(j)j) = 



and U"{(f)i) > ; and one maximum at (j)jjj = (j){rjjj) = fJ,/J\ cr |, with U'{(j)iii) = 
and U"{(j)ni) < {U"{(f)) = |^)- There is also an infiexion point at 0// = (pirji) = 
/ji/Js I 0" I, with U'{(f)ji) > but U"{(f)ii) = 0. Since increases monotonically with 
respect to r, this corresponds to three regions I, II, III in space with rj < rjj < rm. 
Moreover, the monotony of together with the relation 0/// = y^ 0// that links 
(pifiii) to 0(r//) involves (using the solution (^), in the first approximation) : 



riii 



> VS rn (6) 



Let us call 0o, generically, a local extremum of U{(f)) or U'{(f)). In the neigh- 
bour region of space, | — 0o |^ 1, equation (|^) may be solved in the weak field 
approximation by linearizing the function U'{(j)) about 0o. This yields 

"^'^ + -? - U'{<Po) - f/"(0o)(0 - 0o) = 7T^ (7) 



dr r dr Itq r 

• in the first region of space (region I) , 0o = 0/ 

• in region II, 0o = 0// 

• in region III, 0o = 0///. 

Besides, as for the Higgs mechanism of symmetry breaking which too involves a 
scalar field and a quartic self-interaction potential, an analogy can be made with a 
well known phenomenon in solid state physics : the Meissner effect which is a phase 
transition of the second kind, between the superconducting to the normal state. 
Here, the field plays the role of the magnetic fiux and f/"(0) plays the role of the 
difference AT = T^. — T between the temperature, T, of the solid and its critical 



temperature, Tc 



In region I, equ(0) reads 



with solution 



which imphes 



(fcf) 2d(j) 2 / , ,^ 1"s ,gs 

' ' ^' ^•^(l_e-('-'^^)A) (9) 



2ro r 
A^ To ,^ ,^ r 



where A = l//i characterizes the dynamical range of the 0-field in region I. 
Clearly, it is necessary that rj = for the solution (^ be consistent with the 
conditions on (p and -£;. 

Let us notice that, for a sufficiently massive 0-field (m,^ > 10v2 Ti/ctq = 
1.8 10^^^ eV/c^, where r© denotes the radius of the central mass), -£: is smaller 
by a factor 1/100 than the same quantity that would be involved by relation 
(^. This means that the 0-field is expelled out from the central region I, so 
that the Pioneer effect is destroyed here. 

This situation appears analogous to the Meissner effect where the magnetic 
flux is expelled out in the superconducting state (AT > 0). That the 0-fleld 
is expelled from region I (f/"(0) > 0), grants that the orbits of the internal 
planets are not modifled. Its signiflcant action on matter is restricted to regions 
II and III. Figure 1 shows the predicted curve y = ap/a^ versus x = r/X for 
region I. 

In region II and about ru, an approximate solution 0° is obtained by solving 
equation (|^) on account that 0o = (pii ', one flnds : 



r, A U' 



III ^2 



C + j^,r-- + ^^^^r' (11) 
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Figure 1: Predicted curve y = ap/a'^ versus x = r/\ for region I. As one can see, 
the limiting condition y < 0.01 is satisfied for x > lOv^- 



and 



#° 



+ 



A U'i 



>ii 



:i2) 



dr 4ro^ ' ^2 ■ 3 
where A and C are constants of integration. Clearly, the extra potential 
U'{(f)ii)r'^/6 will behave like a positive cosmological constant type term. Let 
us notice that the condition of the weak field approximation, | — (pjj |<^ 1, 
is always satisfied as long as | r — r// |^ tq^/ts in as much as the A term is 
neglected. 

In region II below or beyond rjj, an improved solution = 0° + 50 is obtained 
in the first approximation by adding to the previous solution 0° a correction 
term 6(t>. This involves : 



d^ 



dr 



2 



r dr 



-U" 







(13) 



Now, the "curvature" U"{(f)), and hence its mean value, is positive between 



rj and r// but negative between rjj and rju. In the first approximation, one 
may write U"{(j) < (pn) ^ f/"(0/) = 1/A^ and f/"(0 > (pn) ~ -A;^ with k a 
positive constant. We show in the following that A = ^ defines a wavelength 
related to the part of region II beyond r//. So, replacing U"{(f)°) by the value 
— fc^, equation ([13D becomes in the first order approximation : 



, 2 + -- 7^ + A;'50 = 0. (14) 

dr r ar 

The solution of the above equation is of the form 

r) 

(50 = — sinffcr - $//), (15) 

r 

where _B is a constant of integration and $// is a phase offset. Consequently, 
on account of the solution (pJ] ) and the continuity of at the radius ru, the 
first order solution of equation (^ writes in region II beyond rjj : 



A U\<Pn) 2 ^ B 

— -r 1 r H 

4rn r 6 r 



C + -^r--+ ^;''V ^ + -sinfc(r-rjj), (16) 



d0 r, A S 27rr 27r(r - r,,) . 27r(r - r^j) t/'(0//) 

d;^ = ¥7 + ;^ + H [— ^°^ ( A ^ - "" ^ A )] + ^" 

(17) 

Below r//, the solution is of the form : 



2ro^ r 4ro^ r ' 6 

The continuity of at the radius r// involves : 

A^ r 



-r 



s 



^l-e-"^"/^) (19) 



2ro2 ru 

Further, the above solution involves a critical radius r^ at which the solutions 
of both regions I and II are connected. This critical radius is a solution of the 
following equation : 



r, A U'i 



^11)^2 



C + -^r - - + ^^r' = 0. (20) 

4ro'' r 6 

10 



The constant C is determined from relation ( pOf ) by requiring that r^ and A 
vanish whenever one sets M equal to zero. One finds C = — [/'(0//)/6 Tc^ and 
therefore : 

Tc = \/2Cro (21) 

We will neglect throughout the contribution of cosmological constant type 
terms to the dynamics of the ordinary matter at the scale of the solar system 
since this is known to be very small at present epoch. 

In region III, the solution is of the form : 



^111 + j{l-cos[^ir- mi)]} (22) 



which implies 



d(j) D 2ti 27rr 27r 



where D is a constant of integration and A' = 2tt/J\ U"{(j)iii) \ defines a 
wavelength for the 0-field in region III. Hence, the (p-field would have a 
damped oscillatory behavior in the regions of space where U"{(f)) < 0. 

3 Einstein equations 

3.1 The gravitational field sources 

The metric tensor g^y is solution of the Einstein equations 



1 „ SvrG, 



R,u - -R9,u = -^T^u. (24) 



In the presence of the scalar field (p, its right-hand side 

T,y = r;, + T^^) (25) 

incorporates the energy-momentum tensor of the ordinary matter, T°, and the 
energy-momentum tensor of the 0-field itself. Let us emphasize that the scalar field 

11 



considered in this paper is external to gravity (like the electromagnetic field) but 
obeys the equivalence principle (unlike the electromagnetic field). 

3.2 The weak fields approximation 

Let us denote as g°^^ (resp. g°^") the solution of the Einstein equations for = 
and g^j_i, (resp. g^^) the components of the metric tensor in the presence of the 
0-field (all greek indices run over 0, 1, 2, 3 and x° = ct) ; R^^ denotes the Ricci 
tensor, R = g^'^R^y is the curvature scalar (Einstein's summation convention is 
adopted throughout this paper) and the T^^ are the Christoffel symbols. Hereafter, 
whenever we assume spherical symmetry : x^ = r, x^ = 6, x^ = if (for the sake of 
simplicity, for planar motion </? = | in the following), otherwise the x* * denote the 
Cartesian coordinates {i = 1,2,3). Einstein equations rewrite 

R,^ = ^[(r;. - It°9,^) + (Tj^^ - iT^'^g,^)], (26) 

where T° = g'^^T^p is the trace of T;^ and T(<^) = g'^'^T^JlJ is the trace of T^. In the 
weak field approximation, one gets in particular : Tqo — |T*^'^)(7oo = —>^{U{4>)+j Jdcj)) 
and Tq°q — iT°5foo = \pc^ (weak gravitational field approximation). Furthermore, one 
has in the first approximation 

i?oo = ^V^^oo- (27) 

So, we may write : 

.00 = 1 + 2^^^ (28) 

with 

^'V^ = ^<Ui<P)+Jjd<P) (29) 

V^V/v = 4:7iGp and g^Q = 1 + 2Viy/c^, where p is the density of the ordinary matter. 
Derivating partially equ(^) with respect to (p then comparing with equ(^ yields : 



-86 = ^^^-' + ^"'^' ^2°) 
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on account that the derivative -q^ should be bounded even when extrapolated at 
r = 0. 



4 Equation of motion 

The equation of motion of a test body in the presence of the scalar field writes in 
curved spacetime : 

^ + r>-w^ = -dy + ^w^ (31) 

This means that a force term F^ = mc^[d^(f) — [dip / ds)u^] enters in the right- 
hand side of the equation of motion of a test body of mass m in the presence of the 
(/)-field. The first term of the right-hand side, d^cf), is analogous to the electric part 
of the electromagnetic force whereas the second one —{dip/ds) u^ is analogous to 
the magnetic part. Both terms are necessary to satisfy the unitarity of the velocity 
4-vector {u^u^ = 1, hence Un{u'^'Vu)u^ = 0. Equation ( pT]) may be derived from the 
Lagrangian : 

L = -^e-^{g,,uV + l). (32) 

4.1 Motion in weak fields with low velocity 

In the weak fields and low velocity limit, equation (^) simplifies to 

Now, (?** ~ —1 and Fqq ~ —1/2 (7** dgoo/dx^. Hence, the equation of motion rewrites 
in vectorial notation : 

-^ = -VV'„-/c^V^ + -|- (34) 

where we have set / = ti^ — 1. In next section 4.2, we show that / is positive. 
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4.2 Derivation of the long-range acceleration ap 



The projection of equation ([3^ ) above in plane polar coordinates (r, 9) yields, as- 



suming I dr/dt |<^ c ^/f, the radial component of the acceleration vector, 

^ GM_;^. (35) 



r^ dr 

In the low velocity limit, the tangential component of the acceleration vector, ag, is 
equal to zero (conservation of the angular momentum). Clearly, relation (pSJ) is of 
the form : 

ttr = —{a]\f + ap) (36) 

where ajv = ^^ is the magnitude of the Newtonian radial acceleration and ap is 
the radial acceleration induced by the scalar field, 

ap = fc —. (37) 

dr 

Relation (^) applies to any region of space out of the central mass. Besides, since 

d(f)/dr > 0, / must be positive for ap mimics a missing mass gravitational field (see 

section 7 below). Hence, it follows that the radial acceleration induced by the scalar 

field will be directed towards the central mass as observed for Pioneer 10/11, Ulysses 

and Galileo. 

4.2.1 Region I 

In region I and for r ^ A, the scalar field is expelled out and consequently equation 
(|3^ ) simplifies to : 






— = -{l + f,—)VVM (38) 



or equivalently 



—2 = -G Ur, 39 

dt r^ 



where /o = /(O), Ur = r/r is the radial unitary vector and Mhidden = foi^/ro)'^M 
mimics a hidden mass term (so that the true dynamical mass of the Sun differs from 
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its luminous mass, M©, by the amount /o(A/ro)^MQ). It is worth noticing that the 
kind of missing mass which is invoked here mimics a spherical distribution of dark 
matter located within the Sun rather than a solar halo dark matter. In this respect, 
we will distinguish the hidden mass from dark matter. Throughout, hidden mass 
means extra terms involving the 0-field and that mimic a mass term. Both hidden 
mass and dark matter define the missing mass. 

Furthermore, since the 0-field respects the equivalence principle, equation ( pH] ) in- 
volves a maximum shift Z on the frequency of a photon given by : 

Z = (l + /o4)AV^/cl (40) 

Equation (|iO| ) above allows us as yet to put an upper bound on the possible value 
of /o. Indeed, analysis of the data from the tests of local position invariance ("the 
outcome of any local non-gravitational experiment is independent of where and 
when in the universe it is performed") yields a limit /o < 2 10~^ro^/A^ (see C M. 
Will 13). The local position invariance is one of the three pieces of the equivalence 
principle and, since the 0-field respects the equivalence principle, this is a crucial 
test for this field. As we will see further, the 0-field passes the current tests. Indeed, 
one finds that /o is of the order 10"^ (see subsection 5.1.2). In addition, the study 
of the possible effect of dark matter on the motion of the outer planets involves 
that the missing mass within the Sun is necessarily less than V^~^Mq (see Anderson 
et al. 0). Thence, we may conclude that A ^ tq. Clearly, a 0-field of mass 
m^ > 1.8 10^-'^'^ tV l(? passes all the current tests. 

4.2.2 Region II 

Let us neglect for the moment the contribution of the damped oscillations. We 
will also neglect the 0-term in relation (|30|) so that / ^ /o (i.e., we neglect the 
anharmonic terms). Replacing -^ by the expression (^ obtained for region II, 
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relation (p7|) yields : 



/ 2 

ap = a^(l + 2^) (41) 



where we have set 



and 



r'o = VC^o (42) 



«P = YT^ (43) 



foGM 

turns out to be the asymptotic radial residual acceleration. Besides, combining 
relations ( pT]) and (|4^) above yields 



rc = V2 r'o (44) 

and thence 

ap{r = r,) = 2a^ (45) 

which is also the maximum possible value for ap. Relation ( PD may also be derived 
by requiring the continuity of the derivative d(j)/dr at radius Tc (neglecting the A 
term) . 

5 Interpretation of the data 

It has been questionned why the Pioneer effect has gone undetected in the plan- 
etary orbits of the Earth and Mars. Precisely, the Viking ranging data limit any 
unmodeled radial acceleration acting on Earth and Mars to no more than 0.1 10^® 
cm/s^. Indeed, since the Pioneer effect is expected in region II but not in region I, 
there must be some critical radius Vc which allows one to distinguish between these 
two region of space within the solar system. As region II is defined about the radius 
rii, one may reasonably consider that Tc is of the order r///10 and accordingly the 
anomalous acceleration ap should be negligibly small below the radius rc/2. Indeed, 
our estimate of ru, given in subsection 5.2, is in accordance with the estimate of 
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Tc given in subsection 5.1.1 and both estimations corroborate the fact that the Pi- 
oneer effect is neghgibly small below the asteroid belt. Our scalar field, external 
to gravity but which respects the equivalence principle, provides a solution to both 
the anomalous radial acceleration observed on the spacecraft and the absence of a 
comparable effect on the Earth or Mars. 

Indeed, it is worth noticing that all the spacecraft which undergo the Pioneer effect 
were located at radii well beyond the orbital radius of Mars when the data were 
received from them (the closest spacecraft, Galileo and Ulysses, were in the vicinity 
of Jupiter) . 

Moreover, "no magnitude variation of Op with distance was found, within a sensi- 
tivity of 2 10~^ cm/s^ over a range of 40 to 60 AU" . On account of these facts, we 
conclude that the Pioneer effect is a distance effect and a^ is rather asymptotically 
constant within the regions hitherto crossed by the spacecraft. Above all, the scalar 
field approach leads to the same conclusion. 

5.1 Estimate of ap for Pioneer 10/11 using Ulysses data 

To start with, let us recall that no magnitude variation of ap with distance was 
found , within a sensitivity of 2 10~^ cm/s^ over a range of 40 to 60 AU (the data 
analysis of unmodeled accelerations began when Pioneer 10 was at 20 AU from the 
Sun). Thus we may set ap ~ Op for the Pioneer 10/11. Since we need to be given 
at least one point in the curve ap versus r to be able to determine all the parameters 
needed, our strategy will consist to use a piece of information from the Ulysses data 
(the nearest point) to compute the Pioneer 10/11 data (the farthest points). It is 
worth noticing that the piece of information considered by itself gives no information 
on the magnitude of the long-range acceleration of the spacecraft. It is this feature 
that makes the adopted procedure relevant. To compute Op we need to estimate 
first ro and /q. 
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Figure 2: Predicted curve ap/a'p versus the radius (in AU) from 3.9 AU to 60 AU. 
The curve is asymptotically fiat between 40 AU and 60 AU (the damped oscillations 
expected beyond r// have been neglected). 



5.1.1 estimate of tq and Tc 



As one can see, for r = 2r'o, relation (^TD implies that ap = ^a^. Now, this was 



observed for Ulysses in its Jupiter-perihelion cruise out of the plane of the ecliptic 
(at 5.5 AU). This is also consistent with Galileo data (strongly correlated with the 
solar radiation pressure ; correlation coefficient equal to 0.99) if one adopts for the 
solar radiation pressure (directed away from the Sun) a bias contribution to ap equal 
to (—4 ± 3) 10~^ cm/s^. Hence, we conclude that r'o is approximately equal to half 
of Jupiter's orbital radius, that is r'o ~ 2.75 AU and consequently Vc ~ 3.9 AU on 



account of relation (^4]). Let us assume for the moment ( equal to unity ; this leads 
to conclude that tq ^ 2. To AU. Figure 2 shows the shape predicted for the curve 
ap/a'p versus the radius. The plot starts from the radius r^ = 3.9 AU to the radius 
r = 60 AU. 



5.1.2 estimate of /o and derivation of the magnitude of Op 

2 

In our study on the RC of spiral galaxies, we found that /q is of the order ^^^^^^ where 
Vmax denotes the maximum rotational velocity. This seems to be a general order of 
magnitude for this parameter. So, in what follows, we derive an estimation of /o 

2 

using the relation /o ~ ^^in^^ where Vmax is a maximum velocity to be determined for 
the solar system. In the case of interest in this paper, the 0-field under consideration 
though external to gravity is generated from the Sun (or any other star we would 
have considered). Therefore, it seems natural that Vmax should be a typical velocity 
that is related to the matter components of this star and not a peculiar orbital 
velocity. A suitable value for t>^„^ (see Ciufolini Ijl^), perhaps the best for it involves 
thermodynamics parameters solely, is given by the ratio Pc/pc (assuming the perfect 
gas), where Pc and pc denote respectively the central pressure and mass density of 
the star under consideration. Taking the value of T^ given by solar models (Stix jr5[ , 
Brun et al. [||]), the expression t;^^^ = ^ gives for the Sun : /o = (1.72 ±0.04) 10~^ 
and a'^ = (6.8 ± 0.2) 10~^ cm/s^, in good agreement with the recent results which 
give Op = (7.29 ± 0.17) 10"^ cm/s^ as the most accurate measure of the anomalous 



acceleration of Pioneer 10 (Turyshev et al. 0]). Further, the value computed for 
Op may be corrected to Op = (7.23 ± 0.2) 10~^ cm/s^ by identifying A with tq (see 
subsection 5.2 below) : hence, ( ^ 1.07. 

5.2 Damped oscillations and vanishing of ap 

Figure 1 of the paper of Anderson et al. shows an almost harmonic oscillation of ap 
(nothing is said about this by the authors themselves though) for Pioneer 10 which 
starts at the radius rjj = 56.7 ± 0.8 AU with an amplitude apm of the order ^ 10~^ 
cm s~^ (this is derived by comparison with the uncertainty on ap for Pioneer 10) 
and a wavelength A = 2.7 ±0.2 AU the value of which turns out to be quite identical 
to that of To (let us notice by passing that rjj/X is an integer (= 21)). With these 
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Figure 3: Predicted curve {ap — af^)/a^ versus r (in AU), for region II with param- 
eters ru = 56.7 AU, A = 2.7 AU, C = 1 and B = A/16. 

observational data, relation (^) involves, for C ^ 1 and B = A/16, apm — 0.26 10^^ 
cm s~^ between 56 AU and 60 AU as can be seen in figure 3. 

Furthermore, the calculations carried out in subsection 2 lead to predict the 
decline of ap in the form of damped oscillations beyond r = rui- Hence, since r/// > 
v3 r//, we may confidently expect the decline of ap to occur only beyond r = 96.8 
AU. Let us emphasize that the spatial periodicity A involves a temporal periodicity 
Tp = \/vp, where vp is the speed of the spacecraft. Hence, the periodicity of one 
year found for Pioneer 10 has nothing to do with the orbital periodicity of the Earth. 
Indeed, the coincidence just comes from the fact that vp = 2.66 AU/yr for Pioneer 
10 (at least since 1987). As a consequence, Tp should be greater than one year for 
Pioneer 11 since Pioneer 10 is faster mooving than Pioneer 11. 
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6 Derivation of the steady frequency drift using 
the equivalence principle 

Since, the 0-field obeys the equivalence principle, the steady frequency drift may 
be explained in another way than the Doppler effect thanks to this principle (cf. 
Misner et al [|l^). Actually, the steady frequency drift and the corresponding "clock 
acceleration" —at = —2.8 lO^^^s/s"^ shown by the Compact High Accuracy Satellite 
Motion Progam analysis of Pioneer 10 data may also be interpreted as the analo- 
gous of the gravitational redshift linked to the extra potential term Vp = J apdr 
associated to the scalar field. Indeed, the frequency drift ^^ as well as the clock ac- 
celeration -at = ^^ follow from the relation ^ = ^^('-e)-v>(r) _ _^ jr^+a ^^^ 
where r@ denotes the orbital radius of the Earth and r = r^^ ct (one way, as con- 
sidered by the authors) is the distance of the spacecraft from the Earth. Therefore, 
on account that dr = c dt for the photons (one way), one obtains the observed 
relation ap = atc. In this way, the identity ap = atc seems more natural since, in 
this approach, it is indeed the photons received on Earth from the spacecraft that 
are concerned instead of the spacecraft themselves. 

7 Conclusion 

In this paper, we have presented a possible explanation of the "Pioneer effect", 
without being in conflict with the Viking data or the planetary ephemeris. This is 
based on a possible interaction of the spacecraft with a long-range scalar field, 0, 
which respects the equivalence principle. Like any other form of matter-energy, the 
0-field is a gravitational source through its energy-momentum tensor. Conversely, 
its source is the Newtonian potential of the ordinary matter (in this case, the Sun). 
The calculations were performed in the weak fields approximation, with U a quartic 
self-interaction potential. They gave, near the spacecraft, a residual radial acceler- 
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ation directed towards the Sun, with a magnitude ap asymptotically constant (in 
region II). Both ap and the corresponding clock acceleration, a^, computed from 
our formulas are in fairly good agreement with the observed values. Moreover, a 
scalar field of mass m^ > 1.8 10^^^ eV/c^ will be expelled from region I in a way 
quite analogous to the Meisner effect in a superconducting medium. This limits ap 
to no more than 0.1 10~® cm/s^, from the radius of the Sun to r^ = 3.9 AU. It 
is also found that the ap term should be accompanied with damped oscillations in 
the intermediary region between region II and region III. We also predict beyond 
i^iii = 97 AU (that is, about the year 2009 or so, for Pioneer 10) the vanishing of 
ap in the form of damped oscillations. 

Furthermore, the scalar field theory, as developped in this paper, also gives good fits 
for the rotational curves of spiral galaxies as shown in a previous study. Moreover, 
the same field acts at the cosmo logical scales like a cosmological constant. Prelim- 
inary estimations (Mbelek and Lachieze-Rey, in preparation) from the dynamics of 
the external region of the dwarf galaxy DDO 154, actually the sole galaxy for which 



the edge of the mass distribution has been reached (see Carignan & Purton [|T^) 
led to a value Q^ = 0.43 (Ho/100 km s~^ Mpc~^)~^ in fairly good agreement with 
the value ^a = 0.7 deduced from the Hubble diagram of the high-redshift type la 
supernovae (Perlmutter et al. |jl9[, Schmidt et al. ||2^, Riess et al. ^T|, Garnavich 
et al. [g2[) for Hq about 75 km s^V^pc. 
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